Dark energy generated from a (super) string effective action 
with higher order curvature corrections and a dynamical dilaton 
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We investigate the possibility of a dark energy universe emerging from an action with higher-order 
string loop corrections to Einstein gravity in the presence of a massless dilaton. These curvature 
corrections (up to R 4 order) are different depending upon the type of (super)string model which 
is considered. We find in fact that Type II, heterotic, and bosonic strings respond differently to 
dark energy. A dark energy solution is shown to exist in the case of the bosonic string, while 
the other two theories do not lead to realistic dark energy universes. Detailed analysis of the 
dynamical stability of the de-Sitter solution is presented for the case of a bosonic string. A general 
prescription for the construction of a de-Sitter solution for the low-energy (super)string effective 
action is also indicated. Beyond the low-energy (super)string effective action, when the higher- 
\ curvature correction coefficients depend on the dilaton, the reconstruction of the theory from the 

. universe expansion history is done with a corresponding prescription for the scalar potentials. 

I. INTRODUCTION 
Oh 1 

• Einstein equations in their original form, with an energy-momentum tensor for standard matter on the right 
hand side, cannot account for the observed accelerated expansion of our universe. The late time acceleration of the 
l— ~~ '■ universe, which is directly supported by supernovae observations, and also indirectly, through observations of the 
microwave background, of the large scale structure, weak lensing, and baryon oscillations, poses one of the most 
important challenges to modern cosmology. The standard lore aimed at capturing this important effect is related 
to the introduction of the energy- momentum tensor of an exotic matter with large negative pressure (dark energy) 
7-H , in the Einstein equations. The simplest known example of dark energy (for recent reviews, see [J, HI) is provided 
by the cosmological constant. This does not require any adhoc assumption for its introduction, as is automatically 
present in the Einstein equations, by virtue of the Bianchi identities. The field theoretic understanding of A is far 
I/"") ■ from being satisfactory. Efforts have recently been made to obtain A in the framework of string theory, what leads 
to a complicated landscape of de-Sitter vacua. It is hard to believe that we happen to live in one of the 10 100 or 
more vacua predicted by the theory One might take the simplified view that, like G, the cosmological constant 
A is a fundamental constant of the classical general theory of relativity and that it should be determined from large 
scale observations. It is interesting to remark that the ACDM model is consistent with the observations at present. 
Unfortunately, the non-evolving nature of A leads to a non-acceptable fine-tuning problem. We do not know how the 
present scale of the cosmological constant is related to Planck's or the supersymmetry breaking scale; perhaps, some 
■ deep physics is at play here that escapes our present understanding. 

The fine-tuning problem, associated with A, can be alleviated in scalar field models which do not disturb the 
thermal history of the universe and can successfully mimic A at late times. A variety of scalar fields have been 
investigated to this endp], 0, 0|; some of them are motivated by field/string theory and the others are introduced 
owing to phenomenological considerations. It is quite disappointing that a scalar field description lacks predictive 
power; given a priori a cosmic evolution, one can always construct a field potential that would give rise to it. These 
models should, however, not be written off, and should be judged by the generic features which might arise from them. 
For instance, the tracker models have remarkable features allowing them to alleviate the fine-tuning and coincidence 
problems. Present data are insufficient in order to conclude whether or not the dark energy has dynamics; thus, the 
quest for the metamorphosis of dark energy continues. 

The other alternative for getting accelerated expansion is related to modifications of the geometry itself or the 
left hand side of the Einstein equations. There arc several ways of modifying gravity (for a review, see [6]). Higher 
dimensional (including stringy) effects might lead to large-scale modifications of gravity. Another approach, which is 
largely motivated by phenomenological considerations, is related to the modification of the form of the gravitational 
action (like F(R) gravity, etc). The third intriguing alternative is provided by the higher order curvature corrections 
to Einstein gravity due to low-energy (super)string effective action Q. The leading order correction in the string 
expansion parameter a' is given by a Gauss-Bonnet term which has several remarkable features and which was 
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proposed as a dark energy model [8(. The next-to- leading corrections, cubic and quartic in the curvature, crucially 
depend upon the type of string model under consideration. The higher-order curvature invariants are coupled to scalar 
(dilaton/modulus) fields. One might try to fix these fields by invoking some non-perturbative mechanism. In that case, 
the GB term does not contribute to the four-dimensional equations of motion. However, the higher order curvature 
terms contribute in this case. Their presence crucially modifies the fate of a phantom dark energy universe. Since it is 
difficult to realize any scenario with a fixed dilaton or modulus field, the analysis involving dynamically evolving fields 
becomes very important. A number of papers @, S E3, EH, E3, El, Ejl (see also Ref.Ql]) are devoted to the possibility 
of having dark energy with a GB term and a dilaton/modulus field with non-trivial potential. A steep exponential 
potential exhibits a scaling solution which mimics the background (matter/radiation); the nucleosynthesis constraint 
is satisfied provided the slope of the potential is large (lq|. The scaling solution describes a decelerating universe, ft 
is surprising that the GB term can cause a transition from the matter dominated era to a dark energy universe and 
it can also lead to transient phantom energy, provided the slope of the exponential potential and the dilaton coupling 
to the GB invariant are chosen properly [l6[ , or a more complicated choice of scalar potentials is done [T3| ■ In such 
a scenario with the exponential potential, it is quite difficult to satisfy the nucleosynthesis constraint and, secondly, 
the coupling also becomes very large. Since the introduction of the dilaton potential needs assumptions about some 
non-perturbative mechanisms and the massless dilaton naturally arises in the string loop expansion of the low energy 
effective theory, it is important to explore the possibility of a dark energy solution with a massless dilaton. To this 
effect, the second order curvature correction was considered in Ref. @. This, next-to-leading correction contains a 
higher order Euler density which identically vanishes for space-time dimensions less than six; the other remaining 
term is a curvature invariant of order three. The model can lead to a stable dark energy solution, ft is interesting 
to note that the third order correction in a' crucially depends on the type of string theory model. In this paper we 
incorporate string loop corrections up to order three in a' to the Einstein-Hilbert action with a massless dilaton. We 
investigate the cosmological dynamics of the model and explore whether a particular string type is actually sensitive 
to the existence of dark energy. We also outline a general prescription of the construction of the de-Sitter solution in 
presence of higher order curvature invariants coupled to the dilaton field. 

The paper is organized as follows. In section II, we set up the general evolution equations from string effective 
Lagrangian which incorporates curvature corrections, up to order four in R, coupled to a dynamically evolving massless 
dilaton. In section III, we explore the viability of a dark energy solution for models based upon type II, heterotic 
and bosonic strings in the framework of a perturbative string theoretic set up. Section IV is devoted to the stability 
analysis of the de-Sitter solution in the case of a bosonic string model. In section V, we present a reconstruction 
program for a general action with higher-order curvature invariants coupled to dilaton functions. Section VI outlines 
the construction of the de-Sitter solution in a general case. Section VII presents our conclusions and an outlook. 



II. EVOLUTION EQUATIONS 



The process of compactification of the string theory from higher to four dimensions introduces scalar (mod- 
uli/dilaton) fields which are coupled to curvature invariants. For simplicity, we shall neglect the moduli fields as- 
sociated with the radii of the internal space. In what follows, we consider the low-energy effective string theory action 
0,0 
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where <f> denotes the dilaton field which is related to the string coupling, R is the scalar curvature, L$ denotes the 
scalar field Lagrangian, and C c encodes the string curvature correction term to the Einstein-Hilbert action 7] 
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where a' is the string expansion parameter, £c , C<?\ and describe the leading order (Gauss-Bonnet (GB) term), 

the second order and thi 
have the following form 



the second order and third order curvature corrections, respectively. The terms C^i \ and in the Lagrangian 
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Here 6b, Sr = 0,1 and 

ft 2 = R 2 - iR^RT + R^pR^, (7) 
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The correction terms are different depending on the type of string theory; the dependance is encoded in the curvature 
invariants and in the coefficients (01,02,03) and 6h, 6b' 

• For the type II superstring theory: (a, 02, C3) = (0, 0, 1/8) and 5r = 6b = 

• For the heterotic superstring theory: (ci, 02, 03) = (1/8, 0, 1/8) and 6h = 1, 6b = 

• For the bosonic superstring theory: (ci, 02, 03) = (1/4, 1/48, 1/8) and 5h — 0, 8b = 1 

The higher order curvature corrections look complicated, in general. However, the analysis become tractable in the 
case of a Friedmann-Robertson- Walker universe. In what follows we will specialize to the case of the FRW metric 
with a lapse function N(t), namely, 

3 

ds 2 = -N{t) 2 dt 2 + a{tf ^2{dx 1 ) 2 , (12) 

i=l 

In the FRW background, the leading and next to leading corrections simplify; they depend upon the lapse function 
and its time derivative N. Since only terms linear in TV contribute to the evolution equations, we shall omit the 
higher powers of the time derivative of the lapse function. We then have the following expressions for the curvature 
invariants 
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Here/ = H + H 2 . 

Varying the action (fTJ) with respect to N, N, we find the modified Friedman equation 

3tf 2 = p c + p , (18) 

where 
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and where £i(0), £2(0); £3(0) come from the first, second and third order correction terms, respectively, and can be 
written as 

U{4>) = a' m e 2m ^ ro=l,2,3. (20) 

Let us consider the scalar field equation of motion derived by varying the action {T]) keeping in mind the perturbative 
string theoretic description (V(<p) = 0) 

4> + 3Hj>- - - ^£(3) = 0, (21) 

The evolution equations look complicated, in general (see Eqs. ([55)1 and ([55)1 in the appendix), and the analysis of 
the cosmological dynamics seems to be a difficult task. We, therefore, take a different route in the search of a dark 
energy solution. We shall consider the following simple solution Q and examine its viability in the present case 

H = <j> = (f>ohx^-, 

t t\ 



for ho > 0, and 



H = ^oln^-A (22) 

t s — t t\ 

when ho < with t\ as an undetermined constant. This solution leads to a constant EOS 

^ = - 1 -|^ = - 1 + i ^ 

which corresponds to dark energy (resp. phantom dark energy) for ho > (resp. ho < 0), de-Sitter solution is obtained 
for ho —>■ 00 

We will next analyze in detail whether the evolution equations, (|68|) and (|69|) . exhibit the given solution (|22|) for 
realistic values of the constants /io, ti, and cf>o- By substituting (|22|) into Eqs. (l68|) and (|69j) . we obtain the algebraic 
equations 

-4 + 3/io0o + A(fto)* + h(ho)X 2 + f 3 (h )X 3 = (24) 



and 



^ - 3fc§ + U(h )X + h{ho)X 2 + Mh )X 3 = 0, (25) 



where X = a' /t\ and f' s are given by the following algebraic expressions 

fi(h ) = 5 HB (l2hl-12h 4 o), 
/a(fto) = * fl (2^-6/»S + 6fe§-4/ig), 

/s(/»o) = C(3) (J/$-36/$ + 99^-108^ + 54h$\ 
/ 45 

+(5 H ( - — hi + 90/$ - 144/$ + 108/$ - 54/$ 

+<5 B (-135/$ + 540/$ - 882/$ + 684/$ - 342/$) , 
U(ho) = S HB (-12/$) , 

/sfao) - (-/$ + 3/$ + ^/$) , 



Mho) = C(3) - 7 /$ + 15/$--/$ + 9/$-9/$ 



?/$ + 15/$-f ' 
.),/ ( ^/$ - 15/$ - 15/tg + 36/$ + 9/$ 

— hi - 90h 5 - 75/$ + 198^ + 57/$ J (26) 

where Shb = 0, 1/2, 1 for type II, heterotic and bosonic string, respectively. 

In what follows we would like to analyze the validity of expressions and (PB"|) for realistic values of the constants 
4>o and ti, corresponding to specific values of ho relevant to dark energy observations 
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III. DARK ENERGY SOLUTION 



We shall first examine the existence of dark energy solutions (|22|) in general and then will specialize to particular 
types of string models. We will be interested in finding out whether dark energy can distinguish amongst the string 
types. The case of the bosonic string will be of special interest. 



A. The general case 

We can combine Eqs. ([M)) and (j2"5)) into a single cubic equation as 

A 3 (h )X 3 + A 2 (h )X 2 + AMX - Qh 2 (l - 3h ) 
where the coefficients of X are given by 

A m (ho) =/m(M + 2(l-3/lo)/3+m(M 771=1,2,3. 
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In the case of m < 3, we always have the analytic formulae for the roots, which will be useful for the interpretation 
of the relation and the contribution from each of the correction terms to the solution. 

The positivity of X, the real root of the cubic equation will impose a restriction on the possible values of ho- The 
real solution for (|27p can be obtained from the cubic root formula, as 



X(h ) = si(foo) + s 2 (h Q ) - 
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A 2 (h Q ) 
A 3 (h ) 



(30) 



1. X{ho) = a'/tl > 



The cubic equation has one real root X provided q 3 



> 0. We have checked numerically that the relation 



q 3 + r 2 > is true for all ho in the region 0.8 < \ho\, for all three string types. Eq. (|29|) should then be used in order 
to find the range of ho such that X(ho) > and the corresponding equation of state parameter iu e // be confronted 
with the observations. 



Q c < 1 



Another important constraint on ho is dictated by the fact that 0g < 1, as the dilaton is a real scalar function. The 
cubic equation (I27|) does not involve <po; it enters into the Friedman equation through p c which encodes all higher 
order curvature corrections. Using the Friedman equation (fT5| & (|2"2"]) , we find 
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(31) 



where Q c is the dimensionless density parameter contributed by the correction terms. The constraint f2 c < 1 is 
equivalent to 0§ > 0, otherwise the dilaton would turn complex and this would put a bound on the possible range 
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FIG. 1: Plot of A3 against ho is shown for Type II string model. 
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FIG. 2: A plot of Q. c against ho for the type II string shows that ho > 0, the density parameter £l c < 1. 



of ho which should be combined with the constraint dictated by the positivity of real root of the cubic equation 
(|27p. The range of ho compatible with the two constraints should then be confronted with the observation on the 
equation of state for the dark energy. The recent analysis of the three year WMAP data combined with the supernova 
legacy survey (SNLS) constraints the dark energy equation of state parameter wde- At 68% confidence level, the 
best fit value wde is given by wde = — 1-06^q'q8- If the flat prior is imposed, the parameter is constrained by 
wde = — 0.97^0 09) which translates into a bound on ho, as ho < —11.11 and ho > 6.67. 

We next turn to the individual string models to find out their viability as dark energies, in view of the aforesaid 
constraints. 



B. Type II string 

The case of a Type II string is simplest to investigate. In this case, we have Ai(ho) = ^2(^0) = and the equation 
(1271) reduces to 

A 3 (h )X(h ) 3 -6h 2 (l~3h Q ) = 0, (32) 

which has the following solution 
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In this case the expression of A 3 simplifies to 

A 3 (h ) = f 3 (ho) + 2(1 - 3h )fe(h ), (34) 

with 

h(h ) = C(3) {^K ~ 36A& + 99/ig - 108^ + 54Ag 

/efao) = C(3) (--h 4 + 15/** - -h% + 9h 7 - 9h 8 j . (35) 

Let us first implement the condition X(ho) — ot' /t\ > 0. The sign of A 3 is important for constraining ho using the 
positivity of X. From Eq. (|33|) . we obtain 



A 3 (h ) = C(3) fyftg - 48/i 6 + 81h 7 - Whl + Uhl\ 



We have plotted A 3 in Fig. 1. The plot shows that A3 (ho) > for ho > and A 3 (/i ) < when ho < 0. Using Eq. 
(|33| . we get the possible region that gives X(ho) > as < ho < 1/3 which always yields w e // > 1. Thus, no viable 
solution exists in this case. Therefore, up to 4 th order corrections in R, the Type II superstring model is clearly ruled 
out as dark energy (for the fl c < 1 case, see Fig. 2). 

C. Heterotic string 

In this case A2Q10) = 0, and A\(ho) and A 3 (ho) are given by 

A 3 (h ) = Mho) + 2(l-3h Q )Mh Q ), 

Ai(h ) = /i(M + 2(i-3M/ 4 (M, 



where 



Aa(M = -—hl-84:ht + 27Qhl-252hl-54hl 



+C(3) ( y/io - 48^ + 81h 7 Q - 18h 8 + 54ft^ 
re 1.515/ig- 141.699/ig + 367.367/iS- 273.637/ig + 10.911/ig. (36) 

In this case, we check the consistency of XQiq) and fl c numerically as the formula (|3"3")) does not seem to hold in this 
case. However, as we remark below, A 3 might still be used as a yard stick for the consistency check 

We plot X(ho) in Fig. 3 which shows that X > provided < ho < 23.68. This constraint should be combined 
with Sl c < 1. The plot of £l c in Fig. 4 tells us that either < ho < 5.04 or ho > 23.68. Thus, the allowed range 
for the parameter ho is < ho < 5.04, which corresponds to W e ff > —0.868. Such a value of the equation of state 
parameter is ruled out by recent WMAP 3 and SNLS survey data. However, the combined data (CMB+LSS+SNLS) 
forces the equation of state to vary as —1.001 < wde < —0.875. This results shows that the heterotic string model 
is marginally compatible with dark energy observations. 

D. Bosonic string 

We now turn to the bosonic string, for which A m ^ for m = 1,2,3. All f s contribute to A m in this case. We 
quote below the expression for A 3 . 

A 3 (h ) = -45^ - 492ho + 1530/i£ - 1416/iq - U2h% 
+C(3) (y/io - 4 8^o + - 18^o + 
ps -35.985^ - 549.699^ + 1627.370^ - 1437.640/^ - 277.089/^. (37) 
In order to the check if X(ho) > and £! c < 1, we display our numerical results in Figs. 5 and 6, which show that 
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FIG. 3: Plot of X against ho shows that < ho < 23.68 is the viable range of ho for heterotic string. 
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FIG. 4: Plot of f2 c versus ho for the heterotic string model. The region given by < 1 shows the allowed range of ho, which 
corresponds to 0.8 < ho < 5.04 and ho > 23.68. 



• a. X(h Q ) > for h < -6.189 or h > 

• b. fi c < 1 for ho > 0.5 

Note that (r 2 + q 3 ) > 0, the condition for the existence of one real root of (j27|) . constraints ho to be ho > 0.8. We 
therefore conclude that the allowed range for h is given by h > 0.8 and this corresponds to —1 < w e // < —0.17. 
The requirement for the dilaton to be real clearly excludes the possibility of phantom energy. It is really interesting 
that the bosonic string responds positively to the requirement of dark energy. 

Before moving to the next section, a remark about A3 is in order. In the case of a Type II superstring, Ai and A2 
vanish identically, leading to S2 = 0. The sign of A3 then becomes important for the consistency check on X(ha). In 
the cases of heterotic and bosonic strings this is no longer true. In these cases we have directly checked the positivity 
of X(ho) using numerical treatments. Interestingly enough, we have found numerically that, for a generic range of 
the parameter ho, it turns out that s\ » S2 telling us that (|33|) still holds approximately for numerical values of ho 
which are of interest to us. We then could analyze the bosonic and heterotic models by checking the sign of A3 as 
we did for the case of the Type II string. In fact, our numerical check shows that we reproduce the exact numerical 
results presented here to a good accuracy. 
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FIG. 5: The plot of X versus ho for the bosonic string shows that X(ho) > for ho < —6.189 or ho > 
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FIG. 6: f2 c is plotted here against ho for the case of the bosonic string. 



IV. THE DE-SITTER SOLUTION AND ITS STABILITY 

It is only in the bosonic case that we have the desired solution (normal dark energy). The general expressions 
for the Friedman equation and the scalar field equation of motion, given in the appendix, are difficult to analyze in 
general. However, in the de-Sitter case the equations get simplified and in what follows we will analyze the stability 
of this solution. 

Following Nojiri et. al @] we define two new variables: 



X = 2-, Y = H 2 a'e 24>1 ^ 
H 



(38) 



With these new variables the Friedman equation and the equation of motion for the scalar field (see appendix) can 
be written as 



dX 
dN 



-3X + ^-Y 

00 



-Y 2 



342 



Y 6 



54C(3) y3 



dY 1 + 1 (X 2 



Y l 



»y -'^-^_ 3C( 3)r-(l-« 

2 0o 0o \2 0o 



(39) 
(40) 



These expressions are much simpler, since H is a constant and its derivatives vanish identically. For 0o = 



-0.01 the 
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fixed points (X C ,Y C ) are 



A c — > 


*c — * 


-25.1415 
-2.4860 
2.4860 
29.4780 
-0.0005 ± 0.0060 l 


0.0574 
0.0062 
- 0.0062 
-0.0680 
-0.0072 ± 0.2080 i 



The first critical point in the table above is relevant for us since the second turns out to be unstable. The third and 
the fourth point give a negative value to the string expansion parameter and are, therefore, not relevant. We next 
examine the stability of the solution around the critical point {X c , Y c ) = (—25.1415, 0.0574). The perturbation matrix 
M. has the following form 



M 



-3 ^(12 + 8r c +[1026- 162C(3)]r c 2 ) 

^ - g(l + [114 - 18C(3)]F C ) £ + ^[57 - 9C(3)] - ^(1 + [171 - 27C(3)]T C ) 



The eigen values of the stability matrix M. are: —2161.75 and —1.7903. Therefore, the critical point is a stable 
node. In general fixed points exist for the range \c\>q\ £ (0,0.05882). 




FIG. 7: The phase portrait of cosmological evolution described by |T} in case of bosonic string restricted to de-Sitter for 
4>o = —0.01. Trajectories starting anywhere in the phase space converge at the stable node (-25.1415, 0.0574). 



RECONSTRUCTION OF THE UNIVERSE EXPANSION HISTORY: BEYOND THE LOW-ENERGY 

STRING EFFECTIVE ACTION 



In this section we study a more general gravitational action where the coefficients of the curvature corrections 
depend on the dilaton. In addition, a scalar potential is added. We then consider the reconstruction of such general 
modified gravity from the universe expansion history, following the technique developed in[l8j]. Let us keep £ m {4>) to 

be general functions of the scalar field 
form: 



(3) 

>. For simplicity, however, we neglect C c . Then, the action has the following 



S 



1 



d Xy/—g 



R 
~2 



(41) 
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(31 

Even if one includes C c , we can reconstruct the action, although the final expression becomes rather complicated. 

(3) 

Neglecting C c , the explicit form of the FRW equation and the scalar field equation follows: 
= -3H 2 + -tf + V{4>) - 2A£_ 1 {<j))H 3 ~ 72&(c/>)H (h 2 + H 
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+24; L >|<>)<{ -<>[[{[{- + //)// + 2 [it 1 -:- H) 3 -18H 2 (h 2 + +511" \ . (42) 
= 4> + 3H<p + V'(4>)-2A^((}>) [H 2 +H\ + 24^(0) ( H 6 + (h 2 + h) 3 ) . (43) 



By combining (|42| and (|43| . we obtain 



V(^(t)) = 3iJ 2 - ^0 2 + 24i?a(i)VK(0 + 72£ 2 (<Mt))# (h 2 + H 

+24&{(f>(t)) U(h 2 + h)hH-2(h 2 + h) 3 + 18 (h 2 + hY H' 2 -5H b } , (45) 



-—-!- + i 2 {cp{t)) { -13H 5 ~ 45H 3 H - 3HH 2 + — - 12 [H 2 + H) H } 
a(t) 4 8 I H V / 

+6&((j)){-nH i H -20H 2 H 2 - 4H 3 + (-5H 3 - 7HH - H^j H - (h 2 + H)h\ . (46) 

Hence, if we consider the theory where £,i{4>) and V(<j>) are expressed in terms of two functions, g(t) and f{4>), and 
an arbitrary £2(0) as follows (compare with reconstruction in the less complicated case when only the first order 
correction is present [jj, Il7t [l9l| ) 

V(4>) = 3g'(fm 2 - 2^ + 24 9 '(m))e^U(4>) + 72 ^VW) {gl[mf + gll[m)) 

+246(0) {6 (g'(fm 2 + g"(f(<P))) g'(fm 9 '"(f(4>)) - 2 ( 5 '(/(^)) 2 + 5 "(JW)) 3 
+l8{g'{f{<P)) 2 +g"{f{<P))) 2 g'{fm 2 ~bg\f{ct ) )f) , (48) 

nm \ 9 vm _ _a_ + m^_ wm))S _ ^ mWm)) 



u(4>) 



-zg'um {9" urn? + - 12 G?W)) 2 + sfvm) g"vm} 

+66(0) {-1V(/(0))V(/W) - 2O 5 '(/(0)) 2 .g"(/(0)) 2 - 4g"(/(0)) 3 

+ (-5g'(/(0)) 3 - 7 9 '(/(W(/W) - 5 "'(/W>))) </"(/W) - (s'(/(0)) 2 + <?"(/W>))) <?""(/W)}] ■ (49) 
then it is not hard to check that a solution is given by 

H = g(t), <j> = f~ 1 (t). (50) 

Here / _1 (i) is the inverse function of f (</>). 
An example of this situation is the following: 



and f((f>) to be properly defined, we obtain 



/(/■)= //„/ - //, 1,. ( 1 ) . (r,:i) 



H(t) =H + ^j-. (52) 



12 



When t is small, H (|52|) behaves as that in a universe with a perfect fluid, with W e ff = — 1 + 2/3Hi, and when t 
is large, H behaves as in the de-Sitter space, where H is a constant. Then, if we choose Hi = 2/3, we find that 
before the acceleration epoch, the universe behaves as matter dominated one with io e // — 0. After that, it enters the 
acceleration phase. 
Another example is: 

g(t) = H Q ln-L-H 1 ln( t -^-l) , (53) 



^0 \ ^0 

which gives 



Here Hq, Hi, and to are positive constants. When t is small, H (|54[1 behaves in a way corresponding to the perfect 
fluid case, with w e ff = — 1 + 2/3Hq- Then, if we choose Ho = 2/3, the matter dominated universe occurs. On the 
other hand, when t ~ to is large, H behaves as in the phantom universe with u> e // = —1 — 2/3H± < — 1 and a big 
rip singularity at t = to will appear. The three-year WMAP data are analyzed in Ref. 20], which shows that the 
combined analysis of WMAP with the supernova Legacy survey (SNLS) constrains the dark energy equation of state 
wde pushing it clearly towards the cosmological constant value. The marginalized best fit values of the equation of 
state parameter at 68% confidence level are given by —1.14 < wde < —0.93, which corresponds to Hi > 10.7 as Hi is 
positive. In the case when one takes as a prior that the universe is flat, the combined data gives —1.06 < wde < —0.90, 
which corresponds to Hi > 25.0. Therefore, the possibility that wjje < — 1 has not been excluded. 
Finally, an additional example is the ACDM-type cosmology: 



9(t) = xtt— -r In 
3(1 + w) 



, .'3(1 + 10) , 
!sinh( v ^ ' (t-t Q ) 



a 



2 = 'fpoao 3 ^ ■ (55) 



Here I is the length scale given by the cosmological constant I ~ (l0~ 33 eV) and to is a constant. The time- 
development of the universe given by g(t) (]55[) can be realized in the usual Einstein gravity with a cosmological 
constant A and cold dark matter (CDM), which could be regarded as dust. The corresponding scalar potentials can 
indeed be written explicitly, but they are quite complicated functions. 



VI. CONSTRUCTION OF THE DE-SITTER SOLUTION FOR A GENERAL EFFECTIVE ACTION 



Let us study the possibility of realizing de-Sitter space from the scalar field equation and the Friedman equation. 
The coefficients c\, Ca, C3, and also Sh, 5b in the action (IS])-© depend on what kind of string theory, that is, bosonic 
string, type II superstring theory, or heterotic string theory, we are considering. Furthermore, these coefficients could 
depend on the details of compactification. Moreover, the suitable compactification often induces a scalar potential. 
Hence, here we consider the conditions for the coefficients which allow the de-Sitter space solution. In other words, 
we assume the possibility of a more general effective action like in previous section. 

In the de-Sitter space, the Hubble rate if is a constant 

H = H , (56) 
and all the curvatures are covariantly constant. We also assume the scalar field to be a constant: 

<!> = Po ■ (57) 
For simplicity, C3 terms are neglected and the scalar potential V(cf>) is assumed to be given by 

V(4>) = K,e- 2< ^° . (58) 

Then the scalar equation has the following form: 

= V + 2Acia'x 2 + 96c 2 a' V , x = H e 2po/4 ' . (59) 
On the other hand, the Friedmann equation is reduced to 



= V + 3x- 12c 2 a' 2 x 3 



(60) 
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By eliminating c 2 from (|59|) and (f60|) . one obtains 

Q = 3V& + 8sc + 8cia'a; 2 . (61) 
On the other hand, by eliminating Vb from (|59[) and (|6T))) . we get 

= -1 + 8cia'a; + 36a' 2 c 2 x 2 . (62) 
And by further eliminating the x 2 term from (|6ip and (|62p . we find 



12a'c 2 Vo + c 1 
8a' (cf - 4a'c 2 ) 



x = - T o = ^-T7T2 ■ ( 63 ) 



This expression for a; is not always a solution of the two independent equations (|59j) and (|60|> . or equivalently (|61j) 
and l[6"2"|). The condition for a; in (l6"2"j) to be a solution can be obtained by substituting the expression |B5]) into ([HI] 
(or equivalently 1(62))): 

= 2W a' {c\ - Aa'c 2 ) 2 + 8 (12a'c 2 V + ci) (c? - 4a'c 2 ) + c x (12a'c 2 F + c x ) 2 . (64) 
In the particular case when Vq = 0, Eq. (f6"4")) has the following form: 

9c? = 32a'c 2 , (65) 

and (|6"3")) gives 

• T ° = 8a' (c\ - Aa'c 2 ) = ' (66) 

Therefore, if c\ < 0, there is a possibility that there could occur a de-Sitter space solution. We should note that, even 
if the solution exists, the Hubble rate H itself cannot be determined uniquely. In fact, since 

H = V^e Po/ ^° , (67) 

by choosing po properly, the value of H itself could be arbitrarily changed. Then the value of H could be determined 
by the initial condition. 

Hence, we have presented in the above the condition to be satisfied by the coefficients of our effective action, which 
leads naturally to a de-Sitter universe. In other words, if this condition is fulfilled, the early-time and (or) late-time 
universe can be inflationary (non-singular) due to stringy effects. It goes without saying that, again, the stability of 
such de-Sitter universe should be checked in each case, as it was done for the bosonic string earlier. 



VII. CONCLUSIONS 



In this paper we have considered string loop corrections to the Einstcin-Hilbert action given by ([3]), with a dynam- 
ical dilaton <j). We have explored the cosmological dynamics of the corresponding modified gravity in the framework 
of a low-energy string effective action. For simplicity, we have ignored the contribution of the background (radia- 
tion/matter) energy density. The higher-order string corrections to gravity, specially the third-order correction in 
a' , crucially depend upon the string type. The evolution equations are quite involved and it is difficult to analyze 
them in general. Taking a different route, we have conjectured a particular solution, H = h$/t, ef> = (f>o\a.t/t\ for 
ho > 0[H = ho/(t s —t),<j) = 4>o ln(t s — t)/t\ when ho < 0]. This solution is important from the dark-energy view- 
point; we have carefully checked its viability by enforcing a consistency check on the parameters of the solution. This 
consistency requirement constraints the range of the parameter ho, which defines the effective equation of state. 

The model based upon a Type II string turned out to be the simplest to investigate semi-analytically. The possible 
range of in this case is given by < ho < 1/3 corresponding to an uninteresting equation of state {Weff > 1)- 
Type II is clearly ruled out because the string expansion parameter cannot be negative. Nevertheless, the situation 
might be improved in presence of matter. It would also be interesting to examine the fate of a phantom universe in 
presence of higher curvature corrections with a massless dilaton. 

In the case of the heterotic string, the dark energy solution exists for ho varying as 0.8 < ho < 5.04, which 
corresponds to w e ff > —0.868. WMAP3 data analyzed with the SNLS survey constraints the dark energy equation 
of state as wde = — 0.97^Q'gg at the 68% confidence level, which puts the heterotic string model under pressure. 
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However, the combined (CMB+SNLS+LSS) data forces the dark energy equation of state parameter to vary as 
— 1.001 < wde < —0.875. Thus, the heterotic string is only marginally compatible with observations. 

Cosmological dynamics based upon the bosonic string turn out to be distinguished amongst all possible string types. 
Indeed, the consistency of the model leads to an effective equation of state given by —1 < w e ff < —0.17, which is 
clearly compatible with data. The stability analysis is difficult to carry out in this case, in general. We have studied 
the de-Sitter case (ho — ► oo) in the bosonic case separately and demonstrated its stability. 

A more general action, with higher-curvature corrections coefficients depending on a dilaton was also considered. 
The general reconstruction method could be developed for such theory, so that a realistic universe expansion history 
can be obtained within some class of scalar potentials. An example which proposes a matter-dominance era before 
cosmic acceleration (quintessence, phantom era or ACDAI cosmology) is presented. The de-Sitter universe in such a 
general theory (as well as for the bosonic string) can arise quite naturally. It is known that, with the addition of a 
scalar-Gauss-Bonnct term only, the low-energy string effective action can indeed help in the resolution of the initial 
singularity problem [2l| . The appearance of de-Sitter solution in the general case with higher curvature corrections 
clearly indicates that the resolution of an initial and/or a final singularity of any type (for the classification of future, 
finite-time singularities, see [22]) could be possible taking into account higher-order string loops. 

We conclude that it is not so easy to get the natural dark energy universe from a low-energy string effective action 
(at the very least up to R 4 corrections) . It could turn out that use of even higher order terms is necessary, or that 
the consideration of a different compactification might lead to a more realistic universe. In this respect, it could be 
expected that taking into account stringy non-perturbative effects might help. For instance, one future possibility 
is to consider not only higher curvature corrections but also to include negative powers of such terms (an inverse a' 
expansion?) like in the models with positive and negative powers of the curvature (23j. 
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APPENDIX 

General evolution equations 

The general evolution equations are obtained by varying the action (JTJ) with respect to (j> and the lapse function N. 
The scalar field equation takes the following form 

= 4> + 3H<j) 

+— (-48cic/e 2< ^°# 4 - 192c 2 a'V'^ .H' 6 + c 3 [432((3) - 2736<5 B - 4325 H ] o/V^ H 8 

-48cic/e 2 ^ ff 2 ii - 288c 2 a' 2 e 40/0o i/ 4 ij + c 3 [864C(3) - 5472<5 B - 86AS H ] a' 3 ' ^ H 6 H 

-288c 2 a' 2 e 4 't >/ ' t > H 2 H 2 + c 3 [792((3) - 7056<S B + 1152<f H ] a^e^^H 4 !! 2 - 96c 2 a' V^ff 3 

+c 3 [288C(3) - 4320(5 B - 720<5 H ]a' 3 e 6< ^° H 2 H 3 + c 3 [36C(3) - 1080<S B - 180^]a' 3 e 6<#,/ *°ff 4 | . (68) 
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The Friedman equation in the general case is given by 

3H 2 = \^> 2 + 24c 2 a' 2 e 4< ^ ff 6 - c 3 [72C(3) - 456S B - 72S H ] a' 3 e 6< ^ £T 8 - 432c 2 a' 2 e 40/ *°i/ 4 J ff 

+c 3 [792C(3) - 70565 B - U52S H ] a ,3 e 6< ^ Q ff 6 ij - 288c 2 c/ 2 e 40/ * 1 J ff 2 ij 2 

+c 3 [828C(3) - 10008(5^ - 1656<J ff ] c/ 3 e 6< ^°i? 4 H 2 + 48c 2 a' 2 e 4 ^ ij 3 

+c 3 [168C(3) - 3600<5 S - 600<5s] a' 3 e 6< ^° H 2 H 3 - c 3 [18C(3) - 540<5 S - 90<5s] a' 3 e^^ H 4 

-U4c 2 a' 2 e 4 ' l,/ ' t > H 3 H 3 + c 3 [264C(3) - 2352<5 S - 384<5s] a' 3 1 ^ H b H - 144c 2 a' 2 e 4 ^°HiLff 

+c 3 [288C(3) - 4320<5 S - 720<5s] a' 3 e 6<l>/ ' t > H 3 HH + c 3 [72C(3) - 2160<5s - 360<5s] a' 3 e^^°HH 2 H 

+ — (-48cia'e 2<#,/0o i/ 3 - 288c 2 a' 2 e 4< ^ # 5 + c 3 [864((3) - 5472(5 B - 864(5s] a' 3 e 6< ^°# 7 

-576c 2 a' 2 e 4< ^°# 3 ij + c 3 [1584^(3) - 14112<5 B - 2304<5 ff ] a' 3 ' ^ H 5 H 
-288c 2 a /2 e 4< ^°#ij 2 + c 3 [864C(3) - 12960<5s - 2160(5 ff ] a' 3 1 ^ H 3 H 2 

+c 3 [144C(3) - 4320(5s - 720(5s] a^e 6 ^^ ^^ 3 } . (69) 
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